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1. Part 1. Experimental study of wake-boundary layer 
interaction for different leading edge shapes 

Introduction 

Recent progress in both the linear and nonlinear aspects of stability theory has 
highlighted the importance of receptivity problem [1]. One of the most unclear 
part of receptivity studies is the receptivity of boundary-layer flow to vortical 
disturbances. Some experimental [2] and theoretical [3] results permits to propose 
that quasi-steady outer-flow vortical disturbances may trigger by-pass transition. 
For this reason the experimental and theoretical study of vortex-boundary layer 
interaction is an actual task. In present work such interaction is investigated for 
the vorticity normal to the leading edge. The interest to this type of vortical 
disturbances arises from the theoretical works of M.E. Goldstein et. al. [4-6], 
where it was shown that the flow around the leading edge amplifies them through 
vortex lines stretching. Far downstream, the amplification leads to significant 
distortion of the boundary layer and, finally, to localized flow separation. 

The purpose of experimental part of this work is study of the role of vortex 
lines stretching mechanism in non-uniform flow interaction with boundary layer. 
Such non-uniformity was produced by the wake behind a wire or a set of wires and 
its interaction with boundary layer at the plates of different leading edge shapes 
was investigated. Laminar-turbulent transition caused by wake/boundary layer 
interaction is a topic of investigation too. 

1.1. Experimental setup and equipment 

The experiment was performed in a low-turbulence direct-flow wind tunnel T-36 I 
of Central Aerohydrodynamics institute (TsAGI). The test section is 2.6 rrt long, 
0.5 m wide and 0.35 m high, and is preceded by 12:1 contraction. The free-stream 
turbulence level in the test section is 0.06%, measured in the band 5-1500 Hz at 
velocities greater then 5m/ s. . 

The general outline of experimental setup is shown in Figure 1 a. Interaction 
of wake from vertically stretched wire or a set of wires with the boundary layer 
over the horizontally mounted plate was studied. Four plates with different shapes 
of leading edge shown in Figure 2 were used. These plates of 1S10 mm long, 
500 mm wide and 20 mm thick were made from Plexiglas and had drains for 
static pressure distribution measurements. Plates #1,2 and 4 had semi-elliptical 
leading edges of aspect ratios 8:1, 4:1, and 1:1. The last semicircular leading 
edge is too blunt, so local separated bubble occurred when the flow around it is 
symmetric. In order to achieve the non-separated symmetric flow over possibly 
blunt nose, the leading edge of plate #3 was designed. The contour of this leading 
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edge is described by expression 
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wheie x,y are longitudinal and vertical coordinates of nose contour, H = 
20?7zm - thickness of the plate, a = 4.65. The radii of noses of plates #1,2, 3, 4 are 
1.25mm, 2.5mm, 5.31mm, 10m?n respectively, so the plate #1 had the most 
sharp nose and nose radius grows monotonically with number of plate increase. 

To control the stagnation point position over nose two flaps mounted above 
the plate near its trailing edge were used. The wires were stretched at the holder 
moving horizontally with steps 0.1 mm. All three components of velocity were 
measured with a DISA 55M01 anemometer. A single hot wire probe with hori- 
zontal wire of diameter 5 jum and a sensitive length 0.5 mm was used for stream- 


wise velocity measurements. Spanwise and vertical components of velocity were 
measured by X-wire probe of similar dimensions. The probe was mounted on a 
carriage and was transversed in streamwise and vertical directions. The accuracy 
of probe movement in vertical direction was 0.1 mm. Instead of probe move- 
ment in spanwise dnection the horizontal transmission of the wires was used. 


Coordinate system and general designations used are shown in Figure 1 b. c. 


1.2. Undisturbed flow around plates 

The interaction of non-uniform flow with boundary layer over the plates of dif- 
ferent leading edge shape was studied for free-stream velocities u ^ = 5 m/s and 
17m/ s. the most of results were obtained for Uoo = 17 m/s, so the undisturbed 
flow around the plates for this velocity only will be described in details. Brief 
information concerning flow with = 5 m/s is presented in the end of this 
subsection, for more fully it is described in previous report. 

For plates #1-3 the flow around the leading edge was adjusted to be symmet- 
ric. Symmetric flow over blunt semicircular nose of plate #4 was separated, with 
the local separated bubble located near nose-plate mating appears. To eliminate 
separation the attachment line was moved to upper part of nose and flow around 
plate #4 becomes substantially asymmetric. 

Velocity distributions along the upper and lower sides of nose for all plates 
computed from static pressure measured are shown in Figure 3. For plates #1-3 
the velocity at lower and upper sides coincides within 2% of outer flow velocity, so 
the flow around leading edge was almost perfectly symmetric. Velocity distribu- 
tion over the plate #4 is substantially asymmetric, with velocity at upper side of 
it is smaller then this at the lower side. Moreover, velocities in first two points at 
upper side are approximately the same. So, stagnation point was located between 
two first drains i.e. at y between 2 and 4mm. 

At all plates the flow accelerates initially and then retreads to u near the 
beginning of the plane part of the surface. In accordance with theory, the largest 
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acceleration should occur at the most blunt nose and the velocity maximum 
should decrease as the leading edge becomes sharper. This rule really fulfilled for 
velocity distributions at plates #1-3. The most blunt nose plate #4 stands out 
of this rule, with the velocity maximum at the upper side of it is the smallest. 
The reason of this is the asymmetric flow about leading edge of plate #4. 

Velocity distributions along the plane parts of all plates measured by hot wire 
probe outside the boundary layer (at y = 40mm) are plotted in Figure 4. For 
plates #1-3 there exist a slight negative pressure gradient along x and pressure 
gradient is slightly positive for plate #4. All these gradients are very small, with 
velocity varies by 1% over the length of lm. 

Velocity profiles in the boundary layer at different distances x from the leading 
edge for plates #1-4 are shown in Figures 5-8. Near the leading edge ( x = 20 
and 40 mm) the profiles have maxima at the joint of boundary layer and inviscid 
flow. The velocity at these maxima (shown by open symbols in figure 3) closely 
coincides with velocity obtained from static pressure measurements. Far from 
leading edge velocity profiles tends to Blasius profile for flat-plate boundary layer. 

At plates #1-3 the flow in boundary layer was laminar with r.m.s. velocity 
pulsations in it was less then 0.3% of w,*, over all plates length. Laminar- turbulent 
transition was observed in a boundary layer of plate #4 at x ~ 500?72??r. Two 
factors may cause the transition: positive pressure gradient and amplification of 
outer flow pulsations at the blunt nose. 

Only plates #2 and 4 were tested for outer flow velocity 5 m/s. Velocity 
distributions over noses for this speed were the same as these for = llm/s 
with accuracy ~ 4/%. The distributions of velocity over plane parts of the plates 
for u a o = 17m/s and 5 m/s coincides within accuracy of measurements (~ 0.3%). 
The boundary layer at both plates in flow with u 0 0 = 5 m/s was laminar, with 
r.m.s. velocity pulsations in it were less then 0.2% from u Q 0 . 

1.3. Flow in wakes 

In the reported work four wakes were used as a source of upstream inhomogeneity. 
One of them is steady laminar wake behind a wire of cl = 0.09mm in a flow with 
u oo = 5m/ s. The Reynolds number of this wake based on diameter of the wire 
was Rd — u^djv = 31 and velocity pulsations within it are the same as in the 
free stream. Two next wakes: the wake behind the wire of d = 0.09mm in flow 
with Uco = llm/s and the wake from wire of d = 0.3mm in flow with it, x , = 5 m/s 
were not entirely laminar. Both of them had Reynolds number R d ~ 105, which 
exceeds the critical Reynolds number R » ~ 40 — 50 for steady laminar flow around 
the cylinder [8]. So, the Karman vortex street occurs in the flow behind wire. 
This street decays gradually downstream and for distance greater then 100 d , the 
peak in velocity pulsation spectrum associated with it disappears [8]. Pulsations 
in such decayed Karman vortex streets were about 0.2% of the outer flow velocity. 
It is more then in laminar wake, but it is sufficiently less then in turbulent one. 
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The fourth wake was a wake behind a set of wires of d = 0.09mm placed in 
distance D = 6mm from each other along the span in flow with u <*> = 17 m/s. 

Profiles of mean velocity in all wakes behind single wire were found well 
correlate with formula 


u 






u . . (1.1) 

where u 0 is velocity deficit, b - half-width of the wake, L - distance from the 

wire, A, B, l 0 -constants. For laminar wake it is an exact solution of Prandtl 

equations with A = SaJ&, B = /„ = 0, From experimental data for 

wake behind wire of d = 0.09mm in flow with u TO = 5 m/s it was found that 
A — 1.35. 

For wakes of wire of d = 0.3mm in flow with Uoo = 5 m/s and of wire of 
d - 0.09mm in flow with Uoo = 17 m/s, which are decayed Karrnan vortex streets 
(i.l) is only a fit of experimental data. Velocity deficits it 0 ,and half-widths b of 
these wakes as functions of distance from wire L are plotted in Figures 9 and 
10 respectively. Lines in this figures shows fits (1.1) with constants A, B and l 0 
found from experimental data. 

The flow in the wake from a set of wires is more complicated. Mean velocity 
profiles measured in these wakes for distances from wires L = 60, 170, 270 370 
550, 800mm are shown in Figure 11. Not far from wires ( L = 60m??t) the wakes of 
the individual wires developed independently one from another. As L increases 
( L ~ 1^0, 270 and 370mm) the wakes becomes wider and the interaction of 
neighboring wakes began. Finally, for L = 550 and 800mm the velocity profile 
becomes near-sinusoidal and it remains qualitatively the same for larger distance 
from wires. The flow distortion An, i.e. the difference between maximal and 
minimal velocities in profiles as function of L is shown in Figure 10 together with 
velocity deficit in the wake from one wire. There were observed no distinctive 
difference between the L- dependencies of flow distortion in wakes from the single 
wire and from the set of wires. Theoretically, this difference should exist, but the 
accuracy of the measurements is not enough to observe it. The maximal velocity 
pulsations in the wake from set of wires were about 0.2 - 0.3%, and slightly 
exceeds the pulsations in a single wire wake. 


1.4. Wake-boundary layer interaction 

Interaction of the boundary layer at the plates of different leading edge shapes 
with wakes described in previous section was studied for various distances from 
wire to leading edge. Geometry and flow parameters of all configurations tested, 
that is diameter and amount of wires, velocity, number of plate used, distance 
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from wire to leading edge are listed in the Table. The Table also includes char- 
acteristics of wakes at the leading edge position and other important parameters 
describing these configurations. 

Interaction of wake behind single wire with boundary layer was studied in 
details, so most of this section is devoted to this topic. Only preliminary results 
concerning a wake from a set of wire action on boundary layer are described here. 

1.4.1. Steady boundary layer distortion 

For small outer flow velocity of 5 m/s (configurations 15-19 in Table) the wake 
slightly distorted boundary layer. The maximum velocity deviation from undis- 
turbed flow was only 3 times grater then velocity deficit in oncoming wake. The 
slight effect of vortical disturbances on boundary layer in this regime is caused 
by low unit Reynolds number Ri = u^l/v — 3.45 * 10 s and large boundary- 
layer thickness comparable with the radius of the plate nose. Main results were 
obtained for outer flow velocity 17 m/s corresponding to higher unit Reymolds 
number Ri = 1.17 * 10 6 . For this value of R\ the maximal flow distortion was 
about 10 times grater, then oncoming flow inhomogeneity. 

There were observed three qualitatively different regimes of single wire wake 
/boundary layer interaction: linear regime, symmetric nonlinear regime and anti- 
symmetric nonlinear regime. General features of these regimes will be illustrated 
for interaction of wake behind the wire of d = 0.09mm in flow with u x = 17 m/s 
with plate #2. For this wake and plate combination linear regime was observed 
in configuration 8, when L — 725mm; symmetric nonlinear regime occurred in 
configuration 6 when L = 250mm and antisymmetric nonlinear regime occurred 
for L = 40mm in configuration 4 (see Table). Evolution in x of spanwise distri- 
butions of streamwise velocity in boundary layer observed in these three regimes 
are shown in Figures 12, 13, and 14 respectively. The boundary layer flow dis- 
tortions Aui, as functions of distance from leading edge x for the same regimes 
are plotted in Figure 15. This distortion Aiq, is determined as a difference be- 
tween the maximal and minimal velocities in the spanwise profile measured for y 
corresponding to u — 0.5uoo. 

Linear regime occurs when the distance from wire to leading edge is large 
enough and velocity deficit of oncoming wake is small. In this case the shape 
of spanwise distribution of velocity in boundary layer approximately repeats the 
velocity profile in the wake. The amplitude of boundary layer distortion grows 
almost linearly with distance from leading edge. 

When the wire was placed closer to leading edge, the boundary layer re- 
sponse to the wake becomes noticeable nonlinear. Nonlinearity manifests itself 
in deformation of spanwise profile of velocity in boundary layer (see Figures 13 
and 14). For moderate distance from leading edge L — 250mm, the profile re- 
mains symmetric with respect to the wake centre, so this regime of interaction is 
called as nonlinear symmetric one. If the wire was placed nearby to leading edge 
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( L = 40mm), the boundary layer response to the symmetric wake becomes anti- 
symmetric as shown in Figure 14). This type of wake/boundary layer interaction 
will be called as nonlinear antisymmetric regime. There are two possible anti- 
symmetric regimes: the "right” regime (shown in Figure 14) with maximum of 
velocity at the right side and the ’’left” one with maximum at the left side. Both 
of these regimes were observed and they change each other randomly when the 
wind tunnel was stopped and started again. For definiteness, in all next figures 
dealing with antisymmetric flow distortion the right regime is shown. Unlike the 
linear regime, the flow distortion in both nonlinear regimes reaches maximum at 
about 100mm from leading edge and then decay downstream (see Figure 15). At 
large distance from leading edge (600mm for example), the distortion from weak 
wake initiating linear regime becomes greater than distortion from high-deficit 
wake initiating non-linear regime. 

Streamwise velocity distributions in plane normal to flow direction (y,z— 
plane) for nonlinear antisymmetric, nonlinear symmetric and linear regimes’ are 
shown in Figures 16a), 17a) and 18. In addition to these data, distributions of 
spanwise and vertical components of velocity in this plane for nonlinear anti- 
symmetric and nonlinear symmetric regimes are shown in Figures 16 b) and 17 
b). Vectors plotted in these Figures are proportional to projection of velocity 
onto ( y,z ) plane. Dashed lines y = 2mm shows the outer edge of boundary 
layer where u = 0.99^. For distances from the plate y < 1mm the measure- 
ments of vertical and spanwise velocity components becomes impossible because 
of probe size, so the reconstruction of the flow based on subsequent speculations 
is drawn here. Figure 17 b) shows that nonlinear symmetric regime of boundary 
layer distortion is associated with two counter-rotational vortices located above 
the boundary layer at z = 4mm. In nonlinear antisymmetric regime only one 
streamwise vortex at z = 4mm was observed as shown in Figure 16 b). 

To explain the formation of three different types of boundary layer distortion, 
let s consider the conceptual scheme of wake/leading edge interaction, shown in 
Figure 19. The wake may be considered as a pair of counter-rotational vortex 
sheets. Stretching of these sheets vortex lines around the nose produce a pair of 
counter-rotational streamwise vortices at the plate surface. The lifting of fluid 
between the vortices causes the diminishing of the streamwise velocity, whereas at 
the periphery of them the downwards outer flow leads to growth in the boundary 
lajer velocity. This explains the boundary layer distortion in linear regime. 

In nonlinear symmetric regime the local maximum of velocity in> the middle 
of the wake appears (see Figure 13). Figure 17 a) shows, that this maximum 
becomes more pronounced in the vicinity of the wall. Formation of this maximum 
may be explained if we suppose that when streamwise vortices become strong 
enough two small-size secondary vortices originate within the boundary layer as 
diawn in Figure 1; b). Contrary to main vortices, the secondary vortices drops 
fluid in the middle and produce local velocity maximum. 

The loss of symmetry and antisymmetric regime production may be explained 
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by the instability of streamwise vortex system associated with nonlinear symmet- 
ric regime. This phenomenon is similar to the well-known Crow instability [7] 
of a pair of streamwise counter-rotational vortices which manifests in twisting of 
vortices in a spiral manner. The result ol it should be the lifting ol one vortex 
and the remained vortex will produce an antisymmetric flow distortion. Such 
vortex is really seen in Figure 16 b), where flow in ( y,z )— plane for antisym- 
metric regime is shown. Another point in support of the instability mechanism 
of antisymmetric regime formation provides the flow distortion development in 
configuration 10 shown in Figure 20 a). In this configuration initially symmet- 
ric boundary layer distortion gradually becomes asymmetric with x growth and 
finally reaches completely antisymmetric state at x = 200 mm. 

Contrary to unbounded stream where the Crow instability occurs for all vortex 
strengths, the wall and viscosity should stabilize the vortex pair. The instability 
should occur only if the vortices are strong enough to lift one of them. In accor- 
dance with rapid distortion theory [4], the streamwise vorticity is proportional 
to initial vertical vorticity in the wake which may be estimated as u 0 /b. The in- 
stability should exist if this vorticity is large with respect to some characteristic 
velocity gradient in a flow around the nose. If we suppose the inviscid nature of 
instability, this gradient should be u^/r where r is the radius of nose, and the 
non-dimensional stability criterion I\ = can be constructed. If the insta- 

bility related with vortices/boundary layer interaction is expected, the vorticity 
should be related to velocity gradient in boundary layer over nose. In this case, 

the stability criterion becomes = if- — , where 8 = f— V /2 is the boundary 

I ,1 . t • 1 D “«> \V-CO / * 

layer thickness in the stagnation point. Both criteria for all configurations tested 
are given in Table. Comparison of two criteria shows that transition of a symmet- 
ric flow into antisymmetric some better correlates with inviscid criterion K. As 
a rule, boundary layer flow is symmetric for K < 0.1 and becomes antisymmetric 
for K > 0.2. The exceptions to this rule are configurations configurations 18 and 
19 and 13,14, where the symmetric distortion occurs for K > 0.2. In two first 
configurations (18,19) the unit Reynolds number is too small, so the wake width 
is approximately equal to boundary layer thickness and wake- boundary layer in- 
teraction is fully viscous. In two other configurations 13 and 14 the flow around 
the nose of plate #4 was not symmetric, so the comparison of the results with 
these for symmetric flow around plates #1-3 is not correct. In any case this corre- 
lation shows the significant role of inviscid processes in symmetry /antisymmetry 
transformation. 

Integral characteristics of boundary layer distorted by wake in all three regimes 
of interaction (displacement thickness <5*, momentum-loss thickness 6**, and form- 
parameter H = 8‘/8 m *) as functions of Z are shown in Figure 20 b). This Figure 
shows, that wake strongly deforms the boundary layer and it’s parameters vary by 
an order of unity. In accordance with the previous speculations, in linear regime 
the boundary layer becomes thicker in the wake centre. In nonlinear symmetric 
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regime the displacement thickness grows at sides of the wake and remains almost 
unchanged in it’s centre. The minimum of 8 * in the wake’s centre is probably 
caused by downward flow produced by secondary vortices. The form of 8'(z) de- 
pendance in antisymmetric regime is antisymmetric too, with <5* increases at the 
left side where fluid moves upward and decrease at the right side where downward 
flow ( y,z ) plane was observed (see figure 16 b). 

Vertical distributions of boundary layer distortion A u b measured in different 
configurations and different distances from leading edge x are shown in Figure 
21. For convenience, these distortions are normalized by their maximums and 
y is referred to boundary layer displacement thickness 8 * in the places, where 
the distortions were measured. In. these variables all data are grouped near two 
curves, shown by thick solid and dashed lines, with data for nonlinear symmetric 
and antisymmetric regimes lie near one curve and data for linear regime form 
the second curve. Curve for nonlinear regimes well coincides with profile of low- 
frequency pulsations in boundary layer under enhanced outer flow turbulence level 
measured by Kendall [9]. Maximum of distortion in linear regimes of interaction 
is shifted from the wall to y ~ 1.7(5*. 

To study the role of leading edge in wake/boundary layer interaction the plates 
#1-3 were tested for the same outer flow conditions. Three tests of this type were 
made foi u Q 0 = lTm/s, d — 0.09 mm and distances from wire to leading edge 
L — 40, 250 and 725 mm. Results of them in form of x— dependencies of boundary 
layer distortions A u b and pulsations u at all three plates for L - 40, 250 and 
725mm are shown in Figures 22, 23 and 24 respectively. Distortions A u b plotted 
in parts a) of the Figures were determined as a differences between the maximal 
and minimal velocities in the spanwise profiles measured for y corresponding to 
u = O.Siw Pulsations u shown in parts b) of these Figures are the maximal 
over span values of r.m.s. streamwise velocity pulsations in frequency range 
^ — 1500 Hz measured for the same y. Distortion and pulsations measured in this 
way are close to their maximum values in the section x = const. The data about 
pulsations permits to estimate the influence of pulsations on steady boundary- 
layer distortion development. 

For small ( L = 40mm) and large (L = 725mm) distances from wire to leading 
edge the influence of nose shape on boundary layer distortion is clearly defined. 

In these cases the distortion grows as radius of leading edge increase, with distri- 
butions of distortion over x are qualitatively similar for all shapes of leading edge. 
For L = 725mm (see Figure 24 a)) the linear regime of distortion occurred at all 
plates. As was mentioned above, the distortion for this regime grows with x over 
the entire length of plate. At first glance, it would seem that the distributions of 
distortion over blunt nose plate #3 stands out of this rule, for distortion reaches 
maximum at x = 300mm and then decreases. However, this is not the case for 
the decrease of distortion is caused by laminar-turbulent transition occurred at 
x ~ 300 — 400mm (see Figure 24 b)). Laminar-turbulent transition in this con- 
figuration is considered in subsection 2.4.2 in more detailes (see also Figure 29). 
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At plate #2 the pulsations are sufficiently large to effect on the steady boundary 
layer distortion too. Only at plate the pure steady flow inhomogeneity inter- 
action with boundary layer was observed. For this reason, the data obtained are 
insufficient for finding the quantitative relationship between the radius of nose 
and boundary layer distortion in linear regime. Whether the growth of distortion 
is linear in x or it saturates at some distance from leading edge is unclear yet. 

For all other distances from wire to leading edge ( L = 40 and 250?n?n) the 
boundary layer remained laminar and effect of pulsations on steady distortion 
development was negligible. For L = 40mm nonlinear antisymmetric regime 
occurred at all plates. The distributions of boundary layer distortion in this case 
shown in Figure 22 a) have maxima at moderate distance from leading edge and 
decayed for large x. The maximal over x distortion is approximately proportional 
to radius of leading edge, and distances from leading edge to maximum varies 
in inverse proportion to this radius. It should be remarked that in Goldsteins 
theory [4] the distance from leading edge to singularity is inversely proportional 
to radius of leading edge also. 

Distributions of distortion over x measured for L = 250mm and different 
plates were quite different, because they correspond to different regimes of wake 
/boundary layer interaction. At plate #3 the nonlinear antisymmetric regime 
occurred and distribution of Auj, is similar to that for antisymmetric regime 
shown in Figure 22 a). Pure nonlinear symmetric regime occurred at plate #2 
and distribution of distortion exhibited maximum and then decayed. At the most 
sharp leading edge plate #1 the 2 — profiles of distortion for small x had features 
of nonlinear symmetric regime, but for x > 300mm they becomes quite similar 
to profiles observed in linear regime. So, the distortion initially decays and then 
begins to grow in a manner familiar for linear regime. For small x, as usual, 
the maximal distortion occurred at blunt nose plate ^3. But lor large x the 
distortion at the most sharp nose plate #1 unexpectedly becomes grater, then 
those at more blunt nose plates #2 and #3. 

Boundary layer response on the wake behind a set of 5 wires was studied only 
in four configurations: for L = 40mm at plates #land #3 (configurations 20 and 
23) and for L = 725mm at plates #land #2 (configurations 21 and 22). Re- 
sults in form of x -dependencies of boundary layer distortion A u b and pulsations 
u are shown by dashed lines in Figure 22 (for L = 40mm) and Figure 24 (for 
L = 725mm) together with similar data for wake from single wire. Except config- 
uration 22 with L = 725mm and plate #2, the boundary layer distortion from a 
set of wires differs from those from a single wire less then 15%. In configuration 22 
the distortion from the set of wires is two times greater then this from single wire. 
Consequently, the flow inhomogeneity amplification in this configuration depends 
crucially upon the shape of oncoming velocity profile. Hence, the wake/ bound- 
ary layer interaction in this regime is nonlinear. Figure 25 shows the evolution of 
spanwise profiles of mean velocity and its pulsations measured in configuration 
22. Unexpectedly, no deformation of almost sinusoidal profile of oncoming flow 
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in boundary layer is seen in this Figure. Similar profiles corresponding to inter- 
action of wake from a set of wires with plate #3 for L = 40mm (configuration 
-3) are shown in Figure 26. For such small distance from wires to leading edge 
the oncoming flow is a superposition of five wakes from individual wires. These 
wakes initially interact with leading edge independently, so the velocity profile 
for x - 25mm is a superposition of 5 profiles familiar for single wake /boundary 
layer interaction in nonlinear antisymmetric regime. The orientation of flow dis- 
tortion from these wires is not the same, with three middle wires produced ’’left” 
regime and two side wires excited ’’right” one. The flow distortions in all but one 
wakes are identical and equal to flow distortion from single wire. This explains 
the coincidence of flow distortions. from a single wire and a set of wires for small 
and moderate x. Sufficient interaction between neighboring wakes begins only 
for x > 100mm where laminar-turbulent transition occurred and comparison of 
distortions becomes meaningful. 

Response of boundary layer at sharp nose plate #1 to both near ( L = 40mm) 
and far (L = 725mm) wakes seems to be linear, so the distortions from one wire 
and a set of wires were approximately the same (see Figures 22 and 24). Linear 
behavior of distortion is probably caused by weak vortex lines stretching by flow 
around sharp nose. However, all conclusions concerning boundary layer response 

on wakes linearity drawn here are quite preliminary and should be proved during 
the future work. 

1 . 4 . 2 . Laminar-turbulent transition caused by wake/boundary layer in- 
teraction 

Laminar-turbulent transition excited by wake/boundary layer interaction was 
studied also. Boundary layer distortions A u b and pulsations u as function of x 
measured in configurations where transition takes place are plotted in Figure 27. 
Two different types of transition may be distinguished. The first one is associated 
with origination of high pulsations (~ 6 — 8%) immediately in the first near-nose 
section, where the measurements were made. This type of transition occurred 
in configurations 12, 13, 23 when near wake interacted with blunt nose plates 
#3 and #4. Example of mean velocity and pulsations profiles development for 
this type transition excited by wake from single wire in configuration 12 is shown 
in Figuie 28. Similar results for transition from a set of wires in configuration 
23 may be seen in Figure 26. Possible reason of pulsations origination near the 
leading edge may be local flow separation provoked by wake. Really, the boundary 
layer distortion in the near nose section x = 20 or 25mm in configurations 12, 

23 are very strong (see Figures 28, 26) and the velocity profiles near nose may 
correspond to separated flow in boundary layer. It’s interesting, that in flow 
distorted by a set of wires (configuration 23) the pulsations originated in wake 
behind one of wires (fourth from left side) only. It may be caused by interaction 
of distoitions of different orientations (’’left” and ’’right”) near the nose. Behind 
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the second wire where the same interaction takes place the growth of pulsations 
occurs too, but it begins later at x — 100mm. For regimes of interaction shown in 
Figures 26 and 28 the effect of pulsations on mean flow development is significant 
over whole plate surface. The flow in the centre of wake becomes turbulent at 
small distance ~ 100 — 200mm from leading edge. Further, the turbulent flow 
spreads on the region much wider then the part of boundary layer distorted by 
the wake. 

The second type of transition is characterized by gradual growth of initially 
small pulsations. Complete transition of this type was observed in configura- 
tions 11 and 22 for boundary layer distorted by single wire and a set of wires 
respectively. Development of mean velocity and pulsations profiles in course of 
transitions in these configurations are shown in Figures 29 and 25. For the wake 
from single wire/boundary layer interaction (Figure 29) the maximal pulsations 
initially (at x = 150mm) were observed at the sides of the wake where spanwise 
gradient of velocity is maximal. In this stage of transition pulsations still have no 
effect on mean velocity profile. Further at x = 350mm the additional maximum 
of pulsations appears in the middle of the wake and distribution of pulsations 
becomes wider. Reaching high level of ~ 10% of u <*, pulsations begin to deform 
mean velocity profile with velocity increases at the sides of the wake. Finally, 
at x — 400 mm several maximums of pulsations forming peak-valley structure 
appears and turbulence begins to spread in spanwise direction. In final section 
x — 500mm the turbulent part of boundary layer is about 5 times wider then 
oncoming wake. 

Transition excited by a set of five wires shown in Figure 25 exhibits the same 
features. Initially at x — 300mm maximums of pulsations are located between 
maximums and minimums of mean flow velocity and coincides with maximal 
gradients of velocity in spanwise direction. Further (x = 350mm) maximums of 
pulsations move to the minimums of mean velocity and finally at x = 500mm 
turbulence spreads to width region. It’s interesting, that pulsations at sides of 
disturbed domain are grater then those in it’s central part during entire transition 
process. 

fo reveal the type of disturbances responsible for transition, the power spectra 
of streamwise velocity pulsations were measured. These spectra in configuration 
11 in section x = 150mm, where pulsations began to grow are shown in Figure 
30. From this Figure it is seen that the spectra measured in the centre of wake, 
at side of it, and far from the wake are quite different. A broad wave packet 
with frequencies around 720 Hz can be seen in the spectrum at side of the wake. 
Spectrum in the centre of it contains this packet too, but the major part of its 
energy is confined in low frequency band. The reduced frequency of this wave 
packet F = 10 6 -2tti///u^ o = 227 is approximately equals to the frequency of high- 
frequency pulsations observed in flow with embedded streamwise vortices in [10]. 
The spectrum measured far enough from wake contains packet with f ~ 500 Hz 
or F ~ 150. This packet, in principle, may be formed by Tollmien-Schlichting 
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waves generated in boundary layer by the disturbances within the wake. 

Spectia measured in configuration 22 where transition excited by a set of 
wires was observed are shown in Figures 31 and 32. Figure 31 shows spectra 
measured at the side of wake for z = — 17mm, whereas spectra measured at the 
nearest to centre maximum of mean velocity (z = —3 mm) are plotted in Figure 
32. Three spectra measured at beginning ( x = 150mm), middle (x = 200mm) 
and end (x = 350mm) of transition are presented in both Figures. A packet 
with frequency~ 720 Hz the same as for single wire is seen only in spectrum at 
side of wake in the beginning of transition. Other spectra measured at side of 
wake are broadband, with low frequency pulsations grow predominantly over the 
middle and end of transition. In the centre of wake all spectra are broadband 
and tiansition here is associated with growth of low-frequency pulsations. 

1.4.3. Discussion 

In this section the data obtained will be compared with recent experimental and 
theoretical results dealing with boundary layer receptivity to quasi-steady outer 
flow vortical disturbances and subsequent amplification of these disturbances in 
a boundary layer. For convenience of such comparison, let’s introduce the am- 
plification coefficient k defined as k — Aub/uo. This coefficient as function of 
Reynolds number R = u^x/v computed for different regimes of single wire wake 
interaction with plates #1 and 3 are plotted in Figure 33 a) and b) respectively. 
For the influence of flow inhomogeneity amplitude u 0 is excluded by means of am- 
plification coefficient introduction, the spanwise scale of flow inhomogeneity (6) 
remains as a single parameter describing oncoming flow. From Figure 33 one can 
see that near the leading edge the boundary layer response to flow inhomogeneity 
is almost independent from it’s scale. At large distance from leading edge large 
scale disturbances with b = 1.6mm ( L = 720mm) amplify and small scale ones 
with b = 0.88 and 0.46mm (L = 250 and 40mm) decay or remain almost con- 
stant. Similar conclusion may be drawn from Kendall’s [9, 11] experimental data 
on transition driven by free-stream turbulence if the correlation between the span- 
wise A and time r ~ 1 // scales of free-stream turbulence A ~ const ■ r ~ const / / 
is assumed. Measured in [11] r.m.s. velocity pulsations in six frequency bands 
shown as function of R in Figure 34 reveals the similar behavior as amplifica- 
tion coefficients in Figure 33. Really, near the nose the amplitudes in all bands 
are approximately the same, so disturbances of all scales amplify almost equally 
heie. Foi large Reynolds number (far from leading edge) the low frequency or 
laige scale pulsations grow but small size or high frequency ones remains almost 
independent from R. Another proof of this trend may be seen in Figure 35 where 
broadband r.m.s. pulsations in boundary layer as function of spanwise wavenum- 
* )er ^ = iTTA are plowed. This curve is computed from the cross-correlation data 
given in Figure 5 of [9]. Figure 35 directly shows, that large-size disturbances 
aie moie amplified than the small-size ones. Nevertheless, there exist maximum 
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of amplitude at 0 ~ 2.5 • 10 -4 and very large-size pulsations grow more slowly. 
Estimates of /? in our experiments obtained assuming A ~ 46 are ~ S • I0 -4 
for large scale inhomogeneity corresponding to L = 720 mm and /? = 1.5 ■ 10" 3 
and 3 • 10 3 for near wakes with L = 250 and 40mm respectively. So, the largest 
scale inhomogeneity used in our tests is too small with respect to optimal size of 
disturbances exhibiting maximal amplification in [9]. This fact explains relatively 
weak growth of disturbances observed here with respect to 30 — 60 times ampli- 
fication of low frequency pulsations in boundary layer reported in [11], Another 
possible reason of weak growth of disturbances in our experiment is the influence 
of non-linearity, for boundary layer distortion was very large ~ 0.3 — 0.5?^ with 
respect to r.m.s. pulsation in [11] been 3 - 5%u 00 . 

If we are of the opinion that all scales of flow inhomogeneity used in our 
tests are small and correspond to middle or high frequency pulsations in [9,11], 
the finding of two different vertical profiles of boundary layer distortion shown 
in Figure 21 becomes consistent with data of [11]. It turns out that vertical 
profiles of high frequency pulsations shown in Figure 36 are not self-similar and 
changes with distance from leading edge. For small .r it is similar to well-known 
distribution for low-frequency pulsations shown by (X) in Figure 21 and it s 
maximum moves from wall as x increase. Similar correlation may be seen in 
Figure 21 where dashed curve with near-wall maximum is composed from data 
measured at x = 150mm, and solid curve with maximum far from wall is formed 
by distortions measured at larger x. However, the shift of maximum in linear 
regime may be caused by non-linearity i.e. large amplitude of distortion ~ 30% 
in our experiment. Really, it is the linear regime produces the maximal deviations 
of the boundary layer displacement thickness in the disturbed domain. In the 
centre of it 8 * exceeds the undisturbed value by factor of 1.4 (see Figure 20 b). 
So, it is unclear, what value of displacement thickness should be used for scaling 
of y in this regime. If some averaged 8* ~ 1.28^ (6q is displacement thickness 
of undisturbed boundary layer) was used, the curves for linear and nonlinear 
regimes would coincide. 

The amplification coefficients obtained here may be more directly compared 
with theoretical results of [12], where the steady spanwise-periodic disturbances 
experienced maximal spatial growth were found. The amplification coefficient 
may be directly computed from the results of this work a s k = \/G , where G is 
the energy growth. If we assume that energy in initial section is proportional to 

the amplification coefficient computed in this manner will be proportional to 
that measured here. Figure 37 shows amplification coefficient as function of R 
for various /? computed from results of [12]. This Figure is quite similar to Figure 
33 a) where our experimental results are plotted, with large scale disturbances 
grows with R and small size ones decay in both Figures. Theory [12] gives the 
expression for wave number providing maximal growth as ~ 0.45/ VH. For 
R = 10 6 corresponding to x ~ 800mm in our experiment it gives f3 — 4.5 • 10 -4 
which correspond to spatial period A ~ 13mm or 6 ~ 3.2 mm and is two times 
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grater than maximal spanwise scale used in our tests. 

In general, results dealing with amplification of steady inhomogeneity in 
boundary layer at sharp nose plate #1 obtained here are in qualitative agreement 
with available data about boundary layer response to free-stream turbulence. The 
amplification coefficients at blunt nose plate $3 exhibit the similar properties, 
but they are two times greater then these for sharp nose plate #1. So, the response 
of boundaij' layer at plate ^3 to outer flow turbulence should be two times greater 
then this at sharp nose models usually used in such experiments. Consequently, 
one may suppose, that the similar development of transition at plate #3 will 
take place for twice smaller outer flow turbulence level. Due to most of bodies 
of practical use (wings, turbine blades and others) have blunt leading edges, this 
conclusion may be important for transition prediction. 


2. Part 2. Theoretical study of boundary layer receptiv- 
ity to steady outer flow inhomogeneity and resulting 
laminar- turbulent transition. 

Introduction 

Response of steady boundary layer on the plate to the steady outer flow inho- 
mogeneity (or vorticity normal to leading edge) was studied in works of Goldstein 
et al [7], Here it was found that the flow around the leading edge amplify the vor- 
tical disturbances throw the vortex stretching mechanism. Initially these vortices 
excite distortion of boundary layer flow growing linearly in streamwise direction. 
Far downstream the amplification result in appearance of flow inhomogeneity of 
finite amplitude. In the experiments described in section 1. the finite amplitude 
flow inhomogeneity in boundary layer was observed too. 

Nevertheless, the laminar- turbulent transition is not caused directly by steady 
flow inhomogeneity, but it is provides conditions for flow instability with respect 
to unsteady pulsations. The stability of such boundary layer with steady spanwise 
modulation of velocity profile is studied theoretically in section 2.1. 

The action of vortex stretching mechanism introduced by Goldstein et all 
[6] for the swept wing flow is a topic of section 2.2. In this flow additional 
mechanism of flow inhomogeneity amplification via generation of steady cross- 
flow instability modes was found. Results obtained in this section explains some 
data of expei iment of section 1 . and provides the new area for luture experimental 
work. 

2.1. Stability of boundary layer with steady inhomogeneity of velocity 
profile 

In this section the stability of boundary layer flow with steady spanwise mod- 
ulation of velocity profile is studied theoretically. Such modulation is a model 
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of a streaky structure originating in the boundary layer subjected to free-strearn 
turbulence. The non-uniform flow in boundary layer produced by wake is another 
example of modulated flow. It s well-known that transition is not caused direct] v 
by steady flow distortion, but it s initiated by the growth of high-frequency travel- 
ling waves. For this reason, the stability studies of modulated flows are necessary 
for understanding of transition caused by outer flow turbulence. In contrast to 
generalized Rayleigh s equation approach used in [6], the stability analysis is 
based on more general linearized Navier-Stokes equations. 

2.1.1. Problem formulation 

Consider the stability of boundary layer flow with steady variations of velocity- 
profile in spanwise direction. We’ll use the coordinate sy'stem introduced in sec- 
tion 2 with lengths scaled with averaged over span boundary layer displacement 
thickness 8 . Basic flow Vo(y, r) is considered to be spanwise-periodic with pe- 
riod T 2 and homogeneous in stream wise direction. We present this flow in form 
of Fourier series 


Vo = {f/o(y),0,0}+V t V,- £ V 0n e”^ 

n=-N (2.1) 

Vi = 

where U 0 (y) - Blasius flow profile, V, - flow inhomogeneity, fi = 2n/ r f.. 

There are two general types of unstable disturbances in periodic flows: the 
disturbances of the same period as the basic flow (fundamental disturbances) and 
subharmonic ones. Disturbances of both types will be sought in common form 

e i(ax-ujt) Vp = {u,v,u;} (2.2) 

We’ll consider the temporal stability i.e. complex u> corresponding to real a 
will be sought. 

Substitution of complete flow-field V = V 0 +eV p into Navier-Stokes equations 
and linearization in e gives an eigenvalue problem for u> 

-za)L 2 oV p = L 1 oV f) (2.3) 

where Li and L 2 are linear operators including derivatives with respect to y. 
Discretization of (3.3) in y yields the eigenvalue problem for matrix 

(A— iu>E)F =0 A = D -1 A (2.4) 

__ Here vector F includes discrete representation of disturbances and matrixes 
A and D corresponds to operators Li and L 2 . 


V P = 


N 


E v n (zy 


n——N 
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2.1.2. Numerical method 


In all works where stability of periodic flows is studied the matrixes A and D 
were determined and eigenvalue problem (3.4) is solved. Here we’ll introduce an 
alternative approach to stability of complex geometry flows which is based on 
algorithm for Navier-Stokes equations solution. This approach makes us possible 
to study the stability of any periodic flow using a code for DNS of periodic flows. 
In all of these codes the flow field is presented in form of two-dimensional Fourier 
seiies in (x , z) plane. If we fix a part of harmonics representing basic flow, then 
remaining harmonics will describe the evolution of disturbances. If the amplitude 
of disturbances is small enough, it’s evolution is described by linearized (near the 
basic flow) Navier-Stokes equations. These equations may be written in form 
similar to (3.3) 

^(L 2 o Vj,) = In o Vp 

where Lj and L 2 are the same as in (3.3). All codes for DNS solve the discrete 
version of these equations of form 


^G = AG 

Here A is the same matrix as in (3.4) and G contains discrete representa- 
tion of disturbances. If the finite-difference scheme used in DNS code is known, 
the matrix A may be easily related with matrix of transition B. The former 
matrix describes the evolution of disturbances over one time step r and relates 
disturbances at time t , G k and at time t + r, G as 

G = BG* 


Transition matrix may be directly computed using DNS code. For Crank- 
Nicholson scheme used in our code, matrix B is related with A as 

B = E ~^ A X [ E+ i A 

It is easy to prove that matrixes A and B have the same eigenvectors, and 
eigenvalues of B, y, } are related with eigenvalues of A, \ 3 as 


_ 1 + 

* 1 - Hi 

So eigenvalues uij of stability problem may be obtained from ftj as 

2i 1 — nj 
UJj = - 

T 1 + M j 

The form of disturbances corresponding to in, is determined by eigenvector of 
B associated with pij. 


The method described here was tested by means of computation of Blasius 
flow stability and the secondary instability of finite-amplitude TS wave. The 
results were in excellent agreement with those of traditional methods. 
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2.1.3. Results 


Stability of Blasius boundary layer with R = u^b'/v — 1000 distorted by har- 
monic velocity modulation or by localized flow inhomogeneity were studied. The 
first flow is a model of periodic streaky structure observed under enhanced outer 
flow turbulence level, the second one is a model of the single streak or distortion 
produced by wake studied in section 2. In both flows inhomogeneitv velocity 
vector V, had only streamwise component u t of form 


u i = a f(y)9(z) 


where a is an amplitude, functions f(y) and g(z) defines vertical and spanwise 
distribution of flow inhomogeneity. Vertical distribution for both flows was 


= exp(-^^) - exp(-kffi £) 


max 

y 


exp(— ) _ exp (_l2+wli) 


(2.5) 


with parameters A = 1.2 and y = 1.4 chosen for coincidence with profile of 
low-frequency pulsations measured in boundary layer subjected to free-stream 
turbulence [9]. This profile is plotted by solid line in Figure 22 a together with 
profile of velocity pulsations from [9] shown by points. For the harmonic inlio- 
mogeneity spanwise distribution was 


g(z) — cos /3z 


whereas for localized inhomogeneity it was represented as 


*(*) = 



N 

1 + (1 4- q) ^2 q n cos 

n=l 


q < 1 


( 2 . 6 ) 


The last function with q = 0.63 shown in Figure 22 b has a narrow maximum 
within a = 0 and vanishes over the remanding part of period. 

Due to basic flow symmetry, the disturbances may be symmetric or antisym- 
metric, i.e. the streamwise velocity of disturbances may be even or odd func- 
tion of z. As was mentioned above, the disturbances of both these types may 
have fundamental or subharmonic period. So, there exist four types of unstable 
modes. Figure 23 shows the spanwise distribution of streamwise velocity Rc(u) 
and amplitudes of pulsations |u| for all these modes computed for harmonically 
modulated flow with a = -0.3, a = 0.25, 0 = 0.6. Symmetric and antisymmet- 
ric fundamental modes are shown in Figure 23 a and b, whereas Figure 23 c and 
d shows symmetric and antisymmetric subharmonic disturbances. For symmetric 
modes |u| is largest at the minimums of basic flow velocity u t , whereas maximal 
pulsations of antisymmetric modes coincides with maxima of gradients \dujdz\. 
The phase velocity of symmetric modes coincides with TS wave velocity, but 
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antisymmetric disturbances propagate faster with c = Re(u:)/a = 0.6. Further 
computations reveals that the subharmonic mode is the most unstable one among 
the antisymmetric disturbances and the fundamental mode is the most amplified 
symmetric disturbance. So, only subharmonic antisymmetric and fundamental 
symmetric modes will be considered later. For brevity, these modes will be named 
as antisymmetric and symmetric ones. 

Growth rates 0 = Im(u>) of symmetric and antisymmetric modes as func- 
tions ol amplitude a and spanwise wavenumber (3 are shown in Figures 24 and 
25 respectively. These results were computed for a = 0.25; amplitude depen- 
dance of 6 corresponds to 0 = 0.6, and <?(/?) corresponds to a = 0.3. Growth 
rates ol all modes increase with amplitude growth. For small amplitudes of flow 
inhomogeneity the symmetric modes are the most unstable, whereas antisymmet- 
ric disturbances becomes the most amplified for large amplitudes a > 0.3. The 
(3- dependences of growth rates are different for symmetric and antisymmetric 
modes. Growth rates of symmetric modes are maximal for /3 = 0 and gradually 
decrease with growth of /?. Amplification rates of antisymmetric disturbances 
initially increases with (3 growth, then reaches maximum at (3 = 0.6 and falls off. 
This result contradicts with conclusion of [10] about growth rate of these modes 
proportionality to \du t /dz\. 

Figure 26 shows the growth rates of symmetric and antisymmetric modes as 
functions of reduced frequency F - 10 6 uj/R computed in flow with a = -0.3 and 
(3 = 0.6. Both symmetric and antisymmetric disturbances amplify much rapidly 
then IS waves. The growth rates of these two types of modes are comparable, 
with symmetric modes are the most unstable at large frequency, whereas the anti- 
symmetric modes are the most amplified low-frequency disturbances. Instability 
occurs in wide frequency range 20 < F < 220 with maximal growth rate achieves 
at F = 150. Disturbances in the same frequency range amplified in the boundary 
layer flow with embedded streamwise vortices in experiment [10]. 

Stability localized inhomogeneity flow (3.6) with positive and negative ampli- 
tudes a = ±0.3 was studied in order to reveal the influence of width of inhonio- 
geneity region on flow stability. For this purpose the stability of flow (3.6) with 
q = 0.63 and various (3 with respect to symmetric modes was computed. The 
lesults shown in Figure 2/ demonstrates that the flow with high-speed streak 
has almost the same stability characteristics as homogeneous Blasius flow. In 
flow with low-speed streak the growth rates of disturbances remains the same 
as in Blasius flow if /3 > 0.2. The instability in this flow becomes remarkable 
only if -(3 < 0.2, with growth rate of disturbances increases with (3 diminishing. 

It means that single thin streak does not destabilize flow in spite of inflexible 
velocity profile within it. The remarkable instability in this streak occurs only 
if it's width exceeds a threshold value L* ~ 18*. This fact explains the absence 
of transition in boundary layer distorted by wake studied in section 2. Really, 
in all flow configurations tested, the width of distorted part of boundary layer 
did not exceeds the threshold value predicted by theory. It’s interesting to note 
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that the harmonic modulation of small spanwise period T z < L* (ft > 1) leads 
to remarkable instability (see Figure 25). So, the periodically placed streaks are 
more dangerous that equal sole streak. 

2.2. Receptivity of swept wing boundary layer to small outer flow in- 
homogeneity 

In this section the response of boundary layer on the swept wing to the sinusoidal 
variations of outer flow velocity in spanwise direction is studied theoretically. 
Similar to straight wing, the flow inhomogeneity should amplify in boundary 
layer via Goldstein’s vortex stretching mechanism. The resulting disturbances in 
boundary layer will be quite similar to steady cross-flow instability vortices, so 
these vortices should be generated. This is an additional mechanism of steady 
outer flow disturbances growth in swept wing boundary layer. The simple analytic 
theory describing amplification of outer flow inhomogeneity through the cross-flow 
instability modes generation is developed here. 

2.2.1. Problem formulation 

Consider the viscous incompressible fluid flow over the infinite swept wing. Scheme 
of flow configuration and coordinates systems used are shown in Figure 37. The 
outer flow is assumed to be unidirectional and it’s velocity varies slightly along 
the span. Let’s designate the chord of the wing as X, the mean value of outer 
flow velocity as t/co and the angle between the flow direction and chord as a (see 
Figure 37). The radius of the wing leading edge kL is assumed to be small with 
respect to chord, i.e. k <C 1. However, the wing thickness is not necessary small. 
Further we suppose, that the Reynolds number R = U^L/v is large enough for 
viscous effects been negligeable outside thin boundary layer. To describe flow 
the nondimensional variables are introduced using L and U 0 c as reference length 
and velocity. Two Cartesian coordinates systems shown in Figure 37 (a 1 , y, z) and 
(£,?/, z) will be used further. Origins and applicate axis z of both systems coin- 
cide. Abscisaaxis x of first system coincides with outer flow direction, whereas 
abscisa axis of second one is directed along the chord of the wing. Components 
of vectors in these systems will be designated by subscripts x,y,z and 
respectively. 

Further we suppose that velocity vector V far upstream the leading edge 
(£ — > oo) is presented in form 

V = {V Xi V VJ V z } V x ^l+ eu 00 (y) V y = V z -> 0 £ ^ -oo (2.7) 

where e < 1 is small parameter and function u QO (y) describes the velocity 
variation along the span. At the surface of the wing z = 2 S (£), the no-slip 
boundary conditions should be established 

V(£,7?,z s (0) = 0 (2.8) 
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Boundary conditions (3.1), (3.2) together with Navier-Stokes equations de- 
scribing the flow provides complete problem formulation. In contrast to work of 
Goldstein at al [4], only linear with respect to e solution of this problem will be 
sought here. 

2.2.2. Solution for inviscid flow 

As was mentioned above, the flow around the wing is inviscid outside the thin 
boundary layer whose thickness is ~ R~ 1 / 2 . The solution for inviscid flow is 
sought in form 

V = U 0 + eu 

wheie Uq is a potential basic flow around the wing corresponding to uniform 
outer stream and eu is a disturbance, generated by outer flow inhomogeueity. 
The disturbance is governed by linearized Eulier equations. The solution of this 
equations obtained in [13] can be specialized to present case to give 


u = u* 1 * + V<p 

(2.9) 

u< I )=u TO (y).VA(x,y,x) 

(2.10) 

Aip = -(V, u (I) ) 

(-•11) 

(Vp,n)|„,, KtI|) = 0 

(2.12) 


Here 1 is an integral of equation for streamlines of basic flow 

dx _ dy dz 

Uo X U()y Uq Z 

such, that Y -* y as f — ► — oo, and A (x,y,z) is Laighill’s [2] drift function 


U 0x {x',y(x'),z(x')) J 

Integi ation heie is performed along a streamline of basic flow passing through 
point x, y, z. The difference of drift function between any two points on a stream- 
line is equal to the time it takes a fluid particle to transverse the distance between 
those points. 

Toi subsequent consideration it is more convenient to use coordinates (£, ? /,~) 
fitted to wing chord. In these coordinates the drift function takes form 

A = 7? sin a + A x (£,z) 


A (x,y,z) = x + f 

J — OC 
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where 


Ax(f, z) — £ + J 


cos a 


U, 


- 1 




d( 


(2.13) 


is the drift function corresponding to flow around straight wing with outer 
flow velocity equals 1. To find integral Y , let’s relate dy with d£ tor displacement 
along the streamline. Bearing in mind relationship 

U 0r) , c sino 
dr) = —di = —-dZ 


we find 


dy = cos adrj — sin = sin a 


cos a 


U t 


- 1 




di 


Integration of this along streamline provides the expression for V 


Y = lim y{Z,T){0,z(Z)) = ycosa- A x (£,z)sina 
£—► — 00 


(—-IT) 


Substitution of this into (3.4) and (3.5) gives the expression for u (I) and 
equation for <p in coordinates (£,77,2) 


= u oa (Y) sin a = Uoo (Y)V ± A ± cosa (2.15) 

02 

= — Uoo(y ) • V^A X cos a — u ^ — (V x Ax) 2 ] sin q cos q (2.16) 

where and V x = are projections of and three- 

dimensional operator V onto (£,77) plane. 

Further we introduce the ortohonal coordinate system (.s, //) associated with 
streamlines of potential flow around the wing 


s = 



d$ 

VU^TW^l^const 



d®' 

W 2 + 


where 

= = — = ^ w - — = — = — 

cos a <9£ dz cos a dz DZ 

and T = 0 at the wing surface and $ = 0 in the stagnation point. For 
Jacobian I = ||<?(s, f])/d(Z, J7)|| = 1, the expression (3. 9), equation for <p (3.10) and 
operator V x remain unchanged upon the coordinates transformation. However, 
boundary conditions at the surface (3.6) become more simple 


dip 

dy 


(s,0) =0 


(2.17) 
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Solution for potential p> will be sought in form 

‘r’ — V + <P = cot a J u 00 {Y)dY 

Substitution of this into (3.10), (3.11) gives the following problem for ip 
+ Vi£ = (1 - cot 2 a) u^Y) sin a cos a 

gM)-o (2 ‘ 1S) 

Our task is the finding the solution at the wing surface. To do this let’s find 
the asymptotic behavior of A ± as n -» 0. At first the expression for A ± in 
coordinates ($,^) should be found. For this purpose we relate cM> and <l£ for 
displacement along the streamline as 


= 


d<t> , f d$ , 

+ -Jfedz = Ud£ + Wdz = 


u 2 + w 2 
~u~ d( 


substitution ot this into (3.7) gives 


A±(*,¥) = $ + A 0 (tf)- H f— ll 

v ' Jo U 2 + W 2 J 




where 


W 2 


- 1 


d§ 


_ A » w =/r[z^r 

as shown in [4], A 0 ('F) becomes singular near the wall as 

Ao('I') ~ — rln + Aoo(^) $ 0 

wheie r is radius of leading edge and Aoo('f') remains bounded as 'P — + 0. So, 
the drift function within the wall behaves as 


— > In n + A 0 (s) n -+ 0 


when 


A 0 (s) = $ - r 


1 

L^r 1 


*=0 


d<$> In Uq + const 

a 


Here U 0 (s) is velocity at the wall, a is constant depending from the shape of 
the profile of wing. From this result and (3.8) it follows that 


rt 

rj cos ad — In n sin a — A 0 (s) sin a n — ♦ 0 


(2-19) 


With this asymptotic of Y in hand one can easily show that solution of (3.12) 
within the wall appears in form 


<i P — » n 2 F(s, r),Y ) n — ► 0 
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where F remains bounded when n — > 0. So, V(p vanishes as n — > 0, and the 
solution for velocity of disturbance near the wall in accordance with (3.3), (3.9), 
(3.12) becomes 


(V\ dA± , d<p ' 

■ (y) -aT cos “ + 


Ut) ~ Uoo(Y) sin a + — * (y) as n — > oo 

or) sin a 

u n ^0 

Finally, we have following expressions for components of complete velocity 
vector V near the wall 

V s — ► Uq(s) cos a; V v — > sin a + e U °°^ V n — ► 0 as n — > 0 (2.20) 

sin a 

2.2.3. The boundary layer solution 

Near the wall the boundary layer of thickness ~ R~ 1 ^ 2 should be introduced to 
fit no-slip boundary conditions. In this layer we shall use the rescaled normal to 
wall coordinate Z = ny/R and notations u,v,to for velocity components along 
s .77 and n axis. Due to problem linearity, the boundary layer solution will be 
sought in form 

u = u 0 (s , Z) + e2($, 77, Z) 

u - v 0 (5, Z) + ev(s,7], Z ) 

w = ^S=iuo( 6 , 2) + a2(s,7/, Z) 

P = Po(s,Z) + cp{s,7] y Z) 

where basic flow uo, Vq, w 0 and p 0 correspond to boundary layer in the homo- 
geneous stream and 2 , 2 , 25, p describe perturbations of velocity and pressure 
introduced by outer flow inhomogeneity. Basic flow is described by conventional 
Prandtl equations w T ith boundary conditions 

i£ O (s,0) = uo(«s,0) = ico(s,o) = 0; uq(oc) = Uo(s)cosa t> 0 (oo) = sino 

Disturbances obey the linearized Navier-Stokes equations with boundary con- 
ditions following from (3.13), (3.14) 

2, 2, 25(s, 77, 0) = 0 

u{s,T],z) -» 0; 

as Z — * 00 

u(s, 7 ],z) -► (t) cos q + (Mn Z - A 0 (s) - * In /?) sin a) 


23 



(2.21) 


Foi the sake of simplicity, the sinusoidal flow inhomogeneity of form 

Uoo (y) - 

will be considered further. The period of inhomogeneity is assumed to be of 

the same order as the spanwise period of cross-flow instability vortices so 0 is 
assumed to be of order of 1. ' “ 

To meet the boundary conditions (3.15), (3.16) the solution for disturbances 
ot velocity and pressure will be sought in form 



(s, Z)e ,VH( *< s > + ^) 


( 2 . 22 ) 


where *(s) = -0(A o (s 
follows that amplitudes u, 


Ini?) sin a, 0 = 0cosa. From (3.15), (3.16) it 
v, w should satisfy the boundary conditions 


u{s, 0) = u(s, 0) = w(s, 0) = u(s, oo) = 0 v(s, z) -> — 1 e tfih ' z Z oo 

sin a 

where p = sin a. 

The disturbances (3.17) are the products of amplitudes fi, v, w, p, which vary 
relatively slowly with 5 to fast oscillating exponent. When they are substituted 
into linearized Navier-Stokes equations the terms of order of R, VR and 1 appear. 

the terms of order of 1 were neglected and pressure is eliminated from equations, 
then the following equations for amplitudes of velocity are obtained 

L(a)od) + <5 (M(q)ou) + N(a)o-^ + K(o)oS) = 0; 8 = R~ 1/2 (2.24) 


Here 


tau + i0v + 6— = 0 

os 


L(a) = (zquo -I- i0v o ) ■ 



■_5 2 U 0 . <9 2 u 0 

lOt^^T + 10 


dZ 2 


dZ 2 


s/R\dZ 2 ~ 1 ) 


N(q) -u 0 ^ Z2 - 7 2 j - + 2 ia (iau 0 + i0v o ) (2.25) 

and q = 7 2 = q 2 + 0 2 . Expressions for operators M(a), K(q) describing 

lie effect of flow nonparallelity will not be used further, hence they are omitted 
here. 
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Solution of (3.19) with boundary conditions (3.18) will be sought as a sum 


u 

v 

w 


? 1 

Dr ? (z) + 

W r J 


Go 

G 0 

w 0 


)(s,Z) 


( 2 . 20 ) 


where augend fits the boundary conditions (3.18) and addend satisfy the ho- 
mogeneous boundary conditions 


u o (s,0) = vq(s, 0) — w o (s,0) — G 0 (<s,oo) = Go(^ ? oo) — 0 (2.27] 

One of possible forms of augend in (3.21) is 

V T = [ Z H{Z')e^ lnZ 'dZ' 

sin a sin a Jo 

where H(Z) is an arbitrary function filling the following requirements 
H(Z) — ► 1, Z — y oo; H(Z)/Z 3 — > const Z 0 


This choice permits us to satisfy continuity equation for v rj w r and ensures 
that w ry w r y w r , w r and w* v are finite as Z — > 0. 

Substitution of (3.21) into (3.20) gives the following problem for u 0 , Gu, Wo 

L(a)ow 0 +S(M(a)ow 0 +lS(a)o—^+K(a)ou 0 ) ~ —L(a)ow r — SNl(a)ow r (2.28) 

os 

To solve it, let’s suppose that a is close to eigenvalue of Orr-Sommerfeld 
equation a 0 , corresponding to steady cross-flow instability mode. Physically, it 
means that the direction of cross-flow instability vortex is close to the stream- 
line of undisturbed flow at the outer edge of boundary layer. Introducing small 
parameter 9 we write 


a = cxq + 0A(s ); A(s) ~ 1 

In this case the Orr-Sommerfeld operator L(a) may be presented in form 
L(q) = L(a 0 )+^(a 0 ) ■ A (s) 

and (3.23) becomes 

[L(q„) + ^Ag(ao) + 5 M(q 0 )] ow 0 + <5N(a 0 )o^ + 6K(a 0 )ou 0 + 0(08) + 0(0 2 ) 

— L(a)ottv — 6M(a)ot?7,. 
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If we assume that 0 » 6 the solution of this equation may be sought in the 
frame of multi-scale method as 

wo = ^4(X)$(s, Z) + w x (s,Z) X = js X = 6 - < i (2.29) 

where $($, Z) is an eigenfunction of Orr-Sommerfeld equation corresponding 
to eigenvalue a 0 . Substituting this into (3.22) and retaining only 0(0) terms wo 
obtain 


L(q 0 )oio 1 = -L(a 0 )ouv - AA—(a 0 ) o $ - -^N(a 0 ) o 4> 

This equation for u>i has solution only if right part of it is orthogonal to eigen- 
function <J> + of adjoint to L(q 0 ) operator. The orthogonality condition provides 
the following equation for amplitude A(X) 


dA 

dX 


- iA(XX) ■ A = P( XX) 


where 


P(XX) = - 


($ + , L(Q 0 )ou;y) 
($+,N(a 0 ) o $) 


and ('F,$) _ f 0 fy<PdZ is a definition of scalar product. 

It may be easily shown, that general solution of this equation is 


A = “‘ em + + Of-')! ~ = A(A.V) (2.30) 

where a is constant depending from initial condition. To formulate this con- 
dition, let’s suppose that near the leading edge the difference a - Q(J becomes 
finite, so |A(s)j becomes ~ 1 j6 when s becomes small. Here scaling (3.24) be- 
comes unsuitable and u5 0 ~ 1 or A ~ 0 should be written instead of it. So, if we 
want to match the solution (3.25) with solution near the leading edge we must 
to write 

A ~ 0(9) ~ a -f iP - 0(6) 

Fiom here we have an estimate a ~ 0(9), so first term in (3.25) may be 
neglected and finite expression for amplitude becomes 


A. it 
A 


So, the expression for w in ordinary boundary layer valuables i 

, P(s) 


IS 


W — l- 


a - a 0 


$(s, Z) + 0(a — q 0 ) 
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This expression is valid for R~ l/2 < \a - a 0 | < L If \a - a 0 | becomes 
~ R ~ 1//2 the effect of flow nonparallelism should be taken into account. 

Results obtained reveals that boundary response on outer flow inhomogeneity 
is inversely proportional to difference a - a 0 . To estimate this difference we 
shall use the generalized Faulkner- Scan (F-S) self-similar flow as a simple model 
of boundary layer flow at real swept wing. This flow corresponds to outer flow 
velocity distribution 

Uo = s m cos a V 0 = sin a 

and is defined as 



where functions F and G are solutions of boundary value problem for ordinary 
differential equations 


F" + FF" + (i _ jr'i'j = o G‘" + FG' = 0 

F{ 0) = F'(0) = 0 F'(oo) = 1 G(0) = G'(0) = 0 G'(oo) = 1 

Of course, real flow over the wing should be approximated by F-S Flow with 
exponent m varying with s. Immediately at the nose the flow near stagnation 
point with m = 1 occurs, than flow acceleration becomes weaker and m gradually 
diminishes to zero. After the point of minimum static pressure, the flow begins 
to retard and m becomes negative. So, to estimate the difference a — a 0 over 
all wing surface, it should be computed for F-S flow with various m. Results of 
such computations made for R — 1000, a = 7t/4, s = 1 and several values of 
3 = 0.2, 0.4 ,0.6 are presented in Figure 38. The real part of difference a — Oq 
and it’s absolute value as functions of m are plotted by solid and dashed lines 
in this Figure. For all values of spanwise wavenumber Re(a — q 0 ) is positive for 
positive m and negative if m is negative. As m tends to zero Re(a — a 0 ) goes to 
zero. Due to imaginary part of a 0 being small ~ 10~ 3 and a real, the 

absolute value of a — ao is predominantly defined by it’s real part. So, |a — q 0 | 
becomes small enough as m tends to zero. Unfortunately, cross-flow instability 
modes do not exist for m = 0, so the value of a — a 0 Is not defined in this point 
and minimum of \a — a 0 \ could not be found at all. However, for smallest values 
of m = ±0.001, where computations were performed the absolute values | a — 1 

were about 2 * 10“ 3 , so the boundary layer distortion may be about 500 times 
greater then oncoming flow inhomogeneity. 

In real flow over the wing, the boundary layer velocity profiles with rn ~ 0 
correspond to flow around pressure minimum. So, the maximal boundary layer 
distortion should occur near the minimum of static pressure distribution over 
wing profile. However, this conclusion is quite preliminary due to crudeness of 
F-S flow model. More realistic prediction of boundary layer distortion may be 
reached by means of solution of (3.23) in the boundary layer over real wing using 
PSE- method. 
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3. Conclusions 


i. Large volume of averaged (steady) and pulsational flow characteristics of 
boundary layer //free stream inhomogeneity interaction have been measured with 
plates of different nose shapes and for different free stream inhomogeneities. Flow 

m omogeneities were created by the wakes behind one or several wires stretched 
normal to leading edge. 

2. It was found in a qualitative agreement with Goldstein’s theory [4,61. that 
e deformation of normal to the plate vorticity by the flow around the leading 

edge may be the cause of streamwise vortices (streaks) formation observed in the 
piesence of free-stream turbulence. 

3. T he vortex stretching in course of this deformation leads to significant 
amplification of disturbances in boundary layer. 

4. Three different regimes of single wire wake/ bound ary layer interaction 
re iscoveied. linear, symmetric nonlinear and antisymmetric nonlinear ones 
e regimes are distinguished by different character of downstream amplification 

of disturbances and different spanwise disturbances distribution. The antisym- 
metric regime is the most unexpected as antisymmetric spanwise distribution of 

cmSns 65 m b ° ^ kyer takeS Place for ^mmetric upstream boundary 


5. The increase of wake velocity deficit and decrease of wake width lead 
to transformation of distortion from linear to symmetric nonlinear and then to 

antisymmetric nonlinear regime. The criterion for transition to antisymmetic 
legime was introduced. 

6. Distortion of velocity profile in the boundary layer for nonlinear regimes 

' aVe maximum close to leading edge. This maximum may be the cause of early 
( bypass) laminar-turbulent transition. If not, the distortion then decreases down- 
stream due to viscous dissipation and becomes smaller then gradually increasing 
distortion tor linear regime. 6 

7. The vortex lines stretching and hence disturbance amplification increase 
with increase of plate nose radius. 


8 Interaction of several wakes generated by a set of wires intensify the pro- 
cesses caused by single wake. 

l * 9 / J hC r ake/b ° Undar y layer interaction properties investigated here simu- 
late/reflect boundary layer processes caused by low frequency vorticity of free- 
stream turbulence. 


10. The results obtained permit to suppose that the shape of leading edge 
should have noticeable influence on laminar-turbulent transition caused by outer 
low turbulence/inhomogeneities. Transition Reynolds number for blunt nose 
body may be few times smaller than those observed in well-known experiments 

1 ] , per °™ led W,th shai 'P leadin g edge plates. This difference caused not so 

much by different pressure distribution but by different stretching of vortex lines 
by flow around leading edge. 
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11- The stability of steady spanwise-modulated flow was studied theoretically. 
This flow is a model of streaks found in experiments of Kendall et al [9-11] 
at high tree-stream turbulence level. It was shown, that one isolated streak is 
more stable then periodically placed streaks of the same shape. The isolated 
streak does not destabilize flow at all if it is narrower then threshold value of 
~ 7<5 ‘- Experiment described in Part 1 shows that laminar-turbulent transition 
m periodically modulated boundary layer distorted by the wake from a set of 
wires really occurred earlier then in a boundary layer distorted by single wire 
wake. 

12. The lesponse of boundary layer on the swept wing to the sinusoidal 
variations of outer flow velocity in. spanwise direction was studied theoretically. 
In addition to Goldstein s vortex line stretching mechanism, the amplification of 
outer flow inhomogeneity in boundary layer through the cross-flow instability mode 
generation was found. This phenomenon may lead to extremely high receptivity 
of swept wing boundary layer to outer flow vorticity. 

13. Understanding of streaks formation and their structure opens the hope 
of control the laminar-turbulent transition caused by outer flow vortical distur- 
bances. 
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Table 



* regimes of interaction • 1 - linear, s - symmetric nonlinear, a.s - antisymmetric nonlinear 
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Figure 1. (a) - experimental setup: 1 - plate, 2 - wire(s), 3 - wire holder, 

4 - probe, 5 - wire holder movement, 6,7 - probe movement, 
8,9 - flaps. 

(b) - coordinate system and general designations. 

(c) - parameters of single wire wake: Uq - velocity deficit, 

b - half- width. 
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plate Nfi4 


Figure 2. Contours of leading edges 
of plates Ns 1 - 4 
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Figure 5 . Velocity profiles in boundary layer 
at plate Nsl 
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Figure 6. Velocity profiles in boundary layer at plate Ns 2. 
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Figure 7 . Velocity profiles in boundary layer 
at plate Ns 3 
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Figure 8 . Velocity profiles in boundary layer 
at plate Ns 4 
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Figure 9 . Wake behind wire of d = 0.3 mm 
in flow with = 5 m/s 

a) - velocity deficit u 0 

b) - half width b 
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Figure 10 . Wake behind 1 and 5 wires of d = 0.09 mm 
in flow with u« = 17 m/s 

(a) (•) - velocity deficit in the wake from single wire; 
(*)_- distortion Au/u«> in the wake from 5 wires 

(b) - half-width of wake from single wire b 

























Figure 13. Boundary layer distortion in nonlinear symmetric 
regime of interaction 

(configuration. 17 m/s; plate Nfl2; d = 0.09 mm; L = 250 mm) 










Figure 14. Boundary layer distortion in nonlinear antisymmetric 
regime of interaction 

(configuration 4; u<d=17 m/s; plate Ns2; d=0.09 mm; L = 40 mm) 
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Figure 15. Boundary layer distortion as function of distance from 

leading, edge-far linear (—•—); nonlinear symmetric (— ♦— ) 
and nonlinear antisymmetric (—1—) regimes of interaction. 
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Figurel6. Distribution of three components of velocity vector 
in (y,z) - plane at x = 150 mm measured in 
configuration 4 (nonlinear antisymmetric regime) 

(a) - spanwise profiles of streamwise velocity 

(b) - projections of velocity vectors onto (y, z) - plane 
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Figural 7. Distribution of three components of velocity vector 
in (y,z) - plane at x = 150 mm measured in 
configuration 6 (nonlinear antisymmetric regime) 

(a) - spanwise profiles of streamwise velocity 

(b) - projections of velocity vectors onto (y, z) - plane 
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Figure 18.Streamwise velocity distribution in (y,z) - plane at x=250 mm 
measured in configuration 8 (linear regime) 
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configuration 4; x = 150 mm 
configuration 6; x = 150 mm 
configuration 8; x = 250 mm 
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z, mm 

Figure 20b). Integral characteristics of boundary layer distorted 
by wake as function of z for different regimes 
of interaction. 





+ configuration 3; x = 250 mm 


■ configuration 3; x = 800 mm I 



AUb / A lib max 


Figure 21. Vertical profiles of flow distortion in boundary layer 
measured in different configurations (+■♦•***); 
common curve for linear regimes ( — ); 
common curve for nonlinear regimes ( — ); 
low-frequency pulsations in boundary layer 
subjected to free-stream turbulence [9] (*). 
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Figure 22. Flow distortion AUb (a) and r.m.s. velocity pulsation (b) 
as functions of x for near wakes (L = 40 mm) interaction 
with plates Ns 1 - 3. 

A - plate Nsl ; •- plate Ns2;^ - plate Ns3; 

( ) - single wire wake; (---)- wake from 5 wires. 
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Figure 25 Spanwise profiles of mean velocity u and r.m.s. 

pulsations u' in boundary layer in configuration 22 
(u ro = 17 m/s; 5 wires; plate Ns 2; L = 725 mm) 
















Figure 27. Boundary layer distortion (a) and r.m.s. pulsations (b) 
as functions of x in configurations where laminar- 
turbulent transition occurred. 
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Figure 28 . Spanwise profiles of mean velocity u and r.m.s. 

pulsations u' in boundary layer in configuration 12. 
(Uco = 17 m/s; plate Ns 4; d = 0.09 mm; L= 40 mm) 







Figure 29 . Spanwise profiles of mean velocity u and r.m.s. 

pulsations u' in boundary layer in configuration 11 
(Uoo=17 m/s; plate Ne 3; d = 0.09 mm; L= 725 mm) 
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Figure 31. Evolution of power spectrum streamwise velocity pulsations 
S(f) measured at the side of disturbed area (z = -17 mm) in 
boundary layer distorted by wake from 5 wires in configuration 22 
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Figure 32. Evolution of power spectrum of streamwise velocity pulsations S(f) 
measured in the nearest to centre velocity maximum ( z = - 3 mm) 
in boundary layer distorted by wake from 5 wires in configuration 22 
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Figure 33. Amplification coefficient k= Au b /u 0 as function of Reynolds 

number Re = u^x/v for different wakes interaction with sharp 
nose plate Ne 1 (a) and blunt nose plate Ne 3 (b) 

• - wake from wire of d = 0.09 mm at L = 40 mm, a.s. regime 

♦ - wake from wire of d = 0.09 mm at L = 250 mm, s. and tr. regime 
a - wake from wire of d = 0.09 mm at L = 725 mm, s. linear regime 
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Figure 34. Velocity pulsations measured in different 
frequency renges in [ 9 ] as functions 
of Reynolds number Re = u m x/v. 

1-4 -8 Hz; 2-8 -12 Hz; 3-12 -16 Hz; 
4-16-20 Hz; 5-20 - 24 Hz; 6-24 - 28 Hz 











Figure 36. Vertical profiles of high-frequency (24-28 Hz) velocity 
pulsations in boundary layer measured by Kendall [ 9 ] 
at x= 178 mm (x) and x = 1170 mm (O) 
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Figure 38. 
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Vertical (:,) and ppanwise (b) profiles of flow inho.nogcneily 
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Figure 39. 

Spanwisc distributions of slrcainwise velocity of disturbances lic(U) ( 

and it’s amplitude ( ) in harmonically' modulated flow with 

a=-0.3,d=0.25, ft =0.6 

(a), (b) - symmetric and antisymmetric fundamental modes 
(c), (d) - symmetric and antisymmetric subliarmonic modes 
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Figure 41. Growlli rales of s; 

functions of p in I 



Figure 42. Growth rates v.s. reduced fmniciicj /'• Symmetric ( — + — ) ;u„| 

antisymmetric (— •— ) modes for a=-0J, ( ) . i s , Vi „ 






Figure 43. Growth rales v.s fi for locali/xd inhoiiiogciicit) flow (3.0) with 
( ) - a~0.J; ( ) - (■! ») - \s wave 
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Figure 45. Real part ( ) and absolute value ( ) of 

5 - c*o as functions of exponent m in Faulkner-Scan flow. 

* - 13 = 0.2 +-/? = 0.4 » p = 0.6. 
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